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£>y Abstract 

■ 

The effective average action of Yang-Mills theory is analyzed in the framework of 
exact renormalization group flow equations. Employing the background-field method 
and using a cutoff that is adjusted to the spectral flow, the running of the gauge 
' coupling is obtained on all scales. In four dimensions and for the gauge groups SU(2) 

J> ■ and SU(3), the coupling approaches a fixed point in the infrared. 

o 

■ 1 Introduction and Summary 

Understanding the infrared sector of Yang-Mills theory still represents a challenge in quan- 
tum field theory. The strong coupling of the system and the rich dynamics of its degrees of 
freedom are well beyond the applicability of many field theoretic methods. Even without 
Qh! attempting to solve the theory at one fell swoop, it is already difficult to find (and then 
answer) questions that can be disentangled from the full complexity of the problem. 

In this work, we study Yang-Mills theory in the framework of renormalization group 



X 



(RG) flow equations [I]] for the effective average action 0, concentrating solely on the 
running gauge coupling. Whereas perturbation theory describes asymptotic freedom of 
the coupling in the high-energy limit, it fails to predict anything at low energies except 
for its own failure - manifested by the Landau pole singularity. Even without unveiling 
the complete infrared structure of gauge theories (including confinement and a mass gap), 
an analytic knowledge of the running of the coupling towards lower energies beyond per- 
turbation theory is desirable. Exact RG flow equations represent an appropriate tool for 
tackling this problem. 

Flow equation for the effective average action. Being a "coarse-grained" free-energy 
functional, the effective average action governs the dynamics of a theory at a momentum 
scale k. It comprises the effects of all quantum fluctuations of the dynamical field variables 
with momenta larger than k, whereas fluctuations with momenta smaller than k have not 
(yet) been integrated out. Decreasing k corresponds to integrating out more and more 
momentum shells of the quantum fluctuations. This successive averaging is implemented 
by a fc-dependent infrared cutoff term A^S which is added to the classical action in the 
standard Euclidean functional integral. This term gives a momentum-dependent mass 
square Rk{p 2 ) to the field modes with momentum p which vanishes for p 2 3> k 2 . Regarding 
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Tk as a function of k, the effective average action runs along a RG trajectory in the space 
of all action functionals that interpolates between the classical action S = Tk-^oo an d the 
conventional quantum effective action T = I\.^o- The response of Tk to an infinitesimal 
variation of the scale k is described by a functional differential equation, the flow equation 
(exakt RG equation). In a symbolic notation, 



«9 t r fc = isTr 



d t R k (r® + R k ) 1 , d t = k^-, (i) 



dk J 



(2) 

where T k denotes the second functional derivative of the effective average action with 
respect to the field variables and corresponds to the inverse exact propagator at the scale 
k. 

Flow equation in gauge theories. The use of flow equations in gauge theories, as initi- 
ated in ||, [Q, is complicated by the fact that it is difficult to reconcile the Wilsonian 
idea of integrating out momentum shells of quantum fluctuations with gauge invariance. 
Working with gauge-noninvariant field variables such as gluons and ghosts, a regularization 
of the theory with a momentum cutoff necessarily breaks gauge invariance. Nevertheless, 
gauge-invariant flows can, in principle, be constructed by taking care of constraints im- 
posed by the Ward identities which are modified by the presence of the cutoff Q , , [0 , 
||; in practice, resolving these constraints beyond perturbation theory is highly involved; 
for a review, see ||. 

As an alternative, a formulation in terms of gauge-invariant variables such as, for in- 
stance, Wilson loops may therefore be desirable and so has been proposed and worked 



out in [ 10 1 . Related to this, a gauge-invariant regularization has been formulated in [11 
by constructing SU(iV) Yang-Mills theory from a spontaneously broken SU(iV|iV) super- 
gauge extension; here the fermionic super-partners become massive and act as Pauli-Villars 
regulator fields without breaking the residual SU(iV) gauge invariance. As a result, the 
one-loop (3 function has been computed without any gauge fixing. 

In this work, we decide to employ the conventional and technically more feasible for- 
mulation in terms of the gluonic gauge field at the expense of only partially resolving the 
modified Ward identities resulting in less control over gauge invariance. In this way, we 
shall accept a compromise between calculational advantages and the implementation of 
complete quantum gauge invariance. In particular, we follow the strategy of [0, employ- 
ing the background-field method. Our solution to the flow equation will be gauge invariant 
in the background field, but the renormalization group trajectory that connects the clas- 
sical (bare) action with our quantum solution will not satisfy all requirements of gauge 
invariance (cf. Sect. 0). 

Truncations. Flow equations for interacting quantum field theories can be solved only 
approximately. A consistent and systematic approximation scheme is given by the method 
of truncations. Herein, the infinite space of all possible actions, spanned by the field op- 
erators compatible with the symmetries, is truncated to a subset of operators; the flow 
equation for the complete effective action can then be boiled down to the flow equations 



2 



of the coefficients of these operators (generalized couplings). The renormalization trajec- 
tory in the space of all actions is thereby projected onto the hypersurface spanned by all 
operators of the truncation. For a selected truncation to be able to describe the physics of 
the system, its operators have to cover the dynamics of the relevant degrees of freedom of 
the system under consideration. Since the relevant degrees of freedom in strongly coupled 
quantum field theories such as Yang-Mills theories may change under the renormalization 
flow, a careful and deliberate choice of the truncation is halfway to the solution of the the- 
ory. In view of the many proposals concerning the "true" degrees of freedom in the infrared 
sector of Yang-Mills theory, their systematic study within a flow equation approach would 
be desirable. Along this direction, interesting and promising results have been obtained in 
13] and []14 |, where the choices of the truncation have been based on the monopole picture 



of infrared Yang-Mills theory. 

In the present work, we follow a different strategy: we stick to the "gluonic language" 
and maintain the gauge field as the basic variable. This avoids complications inherent in 



the change of quantum variables, which has to be performed with great care (see, e.g., [15 



and But in order to account for the fact that the "true" infrared degrees of freedom 



may have a complicated gluonic description, we include infinitely many gluonic invariants 
in our truncation; to be explicit, we consider a truncation in which the gauge-invariant 
part of the effective action is an arbitrary function Wk of the square of the field strength 
F, 

r™[ A ] = J w h (0), e := \f« u f° v , (2) 

and the running of the coupling will be extracted from the flow of the linear F® U F® U term in 
Wk, as it is standard in continuum quantum Yang-Mills theory. At weak coupling, it may be 
sufficient to approximate Wk[0] by a finite series, i.e., a polynomial in 9, being justifiable by 
simple power counting (higher operators are suppressed by powers of the ultraviolet cutoff). 
But at strong coupling, those higher operators can acquire large anomalous dimensions that 
completely obstruct a naive power-counting analysis. In fact, our results show that the 
flow of the complete function Wk contributes to the running gauge coupling, and that the 
flow of higher order operators must not be neglected. 

Beyond the approximations involved (i) in choosing Eq. (|2]) as our truncation (and 
neglecting other invariants) and (ii) in resolving the modified Ward identity only partially, 
we make a third approximation (iii) by neglecting any nontrivial running in the ghost and 
gauge-fixing sectors. 

Regulators. For an explicit evaluation of the flow equation, a cutoff function (or regulator) 
Rk has to be specified. This cutoff function is to some extent arbitrary (see Sect. |D|). In the 
denominator of the flow equation (|IJ), it acts as an infrared cutoff for modes with momenta 
smaller than k; its derivative dtRk in the numerator is peaked 5-like around k and thus 
implements the Wilsonian idea of integrating successively over momentum shells. Different 
choices of Rk correspond to different RG trajectories in the space of all action functionals. 
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But by construction, the complete quantum solution T = Tk-^o, being the endpoint of all 
trajectories, is independent of R^. 

This Rk independence of the solution, of course, holds only for exact solutions to the 
flow equation. Approximations such as the choice of a truncation generically introduce a 
cutoff dependence of the final result. On the one hand, this is clearly a disadvantage of 
the method; one is led to study one and the same problem with many different cutoffs in 
order to extract cutoff-independent information. On the other hand, after having accepted 
that exact solutions might never be at our disposal for most quantum field theories, we 
can exploit the cutoff dependence in order to improve our approximations. 

In order to illustrate this point, let us recall that truncations cut a hypersurface out 
of the space of all action functionals. A truncation will be acceptable if the complete 
quantum effective action lies within or close to this hypersurface. But this is not a sufficient 
criterion: imagine a certain exact RG trajectory (corresponding to a certain cutoff function) 
that begins and ends within this hypersurface, but in between develops a large distance 
to the hypersurface. In the exact theory, this flow may largely be driven by operators 
which do not belong to the truncation spanning the hypersurface. Working only within 
the truncation, the contribution of these other operators cannot be accounted for, and the 
so-found solution to the flow will generally be different from the true solution. 

Instead, the optimal strategy would be to choose those exact RG trajectories (and their 
corresponding cutoff functions) that lie completely in (or close to) the hypersurface. But 
strictly speaking, this ideal case is not possible, since the cutoff function generally couples 
the flow to all operators, so that an RG trajectory will never lie only within a restricted 
hypersurface. A more precise criterion would be that the truncated RG trajectory within 
the hypersurface should be equal to (or close to) the exact RG trajectory after projecting 
the latter onto the hypersurface. Then, the flow towards the quantum solution is driven 
mainly by the operators contained in the truncation, and the final result will represent a 
good approximation to the exact one. However, we are currently not aware of any method 
that fully formalizes these ideas. Up to now, the properties of the flow that depend on the 
cutoff function can only be investigated within a given truncation. 

However, a systematic study of cutoff functions has recently been put forward mainly 
within derivative-expansion truncations in scalar and fermionic theories, and "optimized" 
cutoff functions have been proposed [JPT . The optimization criterion focuses on improv- 



ing the convergence of approximate solutions to flow equations; in fact, for scalar O(N) 
symmetric theories, it leads to better results for the critical exponents [18[ . 



Spectrally adjusted cutoff. The class of cutoff functions employed in this work is also 
considered to be improved in the sense mentioned above. In this case, the improvement 
does not refer to the precise shape of the cutoff function, but rather to the choice of 
its argument. Here, we will use not just the spectrum of the Laplace operator (which 
would be the gauge- covariant generalization of the momentum squared), but the full second 

(2) 

functional derivative of the effective average action T k evaluated at the background field. 

The argument of the cutoff function can be understood as a parameter which controls 
the order and size of the momentum shell that is integrated out upon lowering the scale 
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Figure 1: Running coupling a s versus momentum scale k in GeV for gauge group SU(2), us- 
ing the initial value a s (Mz) — 0.117. The solid line represents the result of our calculation 
in comparison with one- loop perturbation theory (dashed line). 

from k to k — Ak. It appears natural that a truncated flow can be controlled better if 
each momentum shell covers an equal part of the spectrum of quantum fluctuations. The 
spectrum itself is not fixed, but k dependent; lower modes get dressed by integrating out 
higher modes. In order to adapt the cutoff function to this spectral flow, we insert the full 

(2) 

T k into its argument, and so obtain a "spectrally adjusted" cutoff. 

This has two technical consequences: first, as the flow equation is evaluated at the 
background field in our truncation, the right-hand side can be transformed into a prop- 
ertime representation; here, we have powerful tools at our disposal that allow us to keep 
track of the full dependence of the flow equation on the field strength squared. Secondly, 
the degree of nonlinearity of the flow equation strongly increases, inhibiting its straight- 
forward analytical or numerical computation even within simple truncations. We solve 
this technical problem by first expanding the flow for the gauge coupling in an asymptotic 
series, and then reconstructing an integral representation for this series by analyzing the 
leading (and subleading) asymptotic growth of the series coefficients. Whereas most parts 
of our work are formulated in d > 2 dimensions and for the gauge group SU(iV), this final 
analysis concentrates on the most interesting cases of d = 4 and N = 2 or N = 3. 

Results. As a result, we find a representation of the j3 function of Yang-Mills theory. For 
weak coupling, we rediscover an accurate perturbative behavior. As the scale k approaches 
the infrared, the coupling grows and finally tends to an infrared stable fixed point, a s —>■ a*. 
Our quantitative results are 

a* ~ 11.3 for SU(2), 

a, ~ 7.7 ±2 forSU(3). (3) 

The uncertainty in the SU(3) case arises from an unresolved color structure in our calcu- 
lation (cf. App. [E|. 
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The complete flow of the running coupling is depicted in Fig. [I] for pure SU(2) Yang- 
Mills theory in comparison with perturbation theory. For illustrative purposes, we use 
a s (Mz) — 0.117 as initial value (Mz — 91.2 GeV). Sizeable deviations from perturbation 
theory occur for k < 1 GeV, and the fixed point plateau is reached for k = O(lOMeV). 
We shall argue below that a larger truncation as well as the inclusion of dynamical quarks 
are expected to decrease the value of a*. 

The paper is organized as follows: Sect. |^ briefly recalls the framework of flow equations 
in gauge theories with the background-field method and describes our basic approximations. 
In Sect. H, we boil down the flow equation as required for our truncation. Sect. £| is devoted 
to extracting the RG flow of the running gauge coupling, which is the main result of the 
present work. The role of the spectrally adjusted cutoff is illustrated in Sect. |[ Sect. || 
contains our conclusions and a discussion of our results in the light of related literature. 



2 Flow equation for Yang-Mills theory 

We begin with a brief outline of the flow equation and the background-field formalism as 
they are employed in this work. We focus on direct applicability and the required approxi- 
mations and leave aside more formal (though important) aspects, as they are presented in 
T2fl and Let us therefore start with a more explicit representation of the flow equation 



for the effective average action, 



d t T k [A, A] = ~ STr { d t R k (I k 2) [A, A}) 



rf\A,A] + R k (r^[A,A}) 



(2)r 



(4) 



where we denote the so-called classical gauge field by A a , which is the usual field variable 
of the quantum effective action (conjugate to the source). We also introduce a background 
field A a , and have already inserted T k evaluated at the background field into the cutoff 
function.[] The symbol STr implies tracing over all internal indices and provides for a minus 
sign in the ghost sector. We aim at solving Eq. (W), using the following truncation: 



r k [A,A] 



Tf[A] + If [A, A] + If [A, A] + rr se [A A]. 



(5) 



Following ||2Q|| , the background-field method is introduced to enable us not only to perform a 
meaningful integration over gauge-fixed quantum fluctuations but to simultaneously arrive 
at a gauge-invariant effective action. Identifying the quantum fluctuations with A — A, the 
gauge-fixing term 



Tf[A,A] = ±J [D ll [A\ (A-A),}'- 



(6) 



1 This Tr' is evaluated at the background field because an A dependence would spoil the one-to-one 
correspondence of the flow equation to the functional integral. 
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with a gauge parameter a is invariant under a simultaneous gauge transformation of A a 
and A a i and so is the ghost action 

TfM = - f ZD»[A]D»[A]c, (7) 

J X 

where the ghosts c, c are understood to transform homogeneously. We should stress that 
with this truncation of the ghost and gauge-fixing sector, we neglect any running there. 

If we solved the theory completely, the resulting quantum effective action r fc=0 [A,A] 
would be gauge invariant precisely at A = A. Imposing a normalization of r| auge [A, A] 
such that 

Tf use [A = A,A] = 0, (8) 

we conclude that the so- found solution r™ L4] would be gauge invariant and would repre- 
sent the desired quantum effective action (provided that we also worked out the complete 
ghost sector). The quantity r| auge hence parametrizes the gauge- noninvariant remainder 
of the action in the physically irrelevant case of A ^ A. 

Although we are finally interested in T k [A = A, A], r| auge cannot be dropped right from 
the beginning in Eq. (f|), because its second functional derivative (r| augc )( 2 ) [A, A] ^ in 
general; rf, auge contributes to the flow of T k [A, A] even at A = A. Still, neglecting r| auge 
seems to be a consistent truncation if we are interested only in P^fA = A]. But besides 
Eq. OH), the effective action also has to satisfy the constraints imposed by gauge invariance 
in the form of the modified Ward identity; in symbolic notation, 

C w [T k ]= A[R k ], (9) 

where £\y denotes the usual Ward operator constraining the effective action 1^, and the 
cutoff-dependent right-hand side represents the modification due to the infrared regulator 
R k (for an explicit representation of Eq. (§), see [0, In fact, if the cutoff is removed 

in the limit k — > 0, we rediscover the standard Ward identity £w[r*;] = 0. Inserting our 
truncation ([5]) into the Ward identity, the first three terms drop out and we are left with 

£ w [r gaugc ] = A[R k }. (10) 

This tells us that, on the one hand, T k nv is indeed not constrained by the modified Ward 
identity and any gauge-invariant ansatz is allowed; on the other hand, a vanishing r| aug0 is 
generally inconsistent with the constraint. It is a nontrivial assumption of this work that 
■pgauge ag ^ r j ven ky right-hand side of Eq. (|10"D does not strongly influence the flow of 
r^ v at A = A, so that we can safely neglect it in a first approximation. 

With regard to our final asymptotic analysis of the running coupling, we can even 
weaken this assumption a bit: since we reconstruct the (3 function from its asymptotic series 
expansion by analyzing its leading growth, neglecting r| auge corresponds to assuming that 
■pgauge ^ Qes nQ ^ gangly modify this leading growth. In view of the fact that r gauge starts 
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from zero in the ultraviolet and enters the flow only indirectly, this assumption appears 
rather natural, at least for a large part of the flow. 

We should remark that the effective average action has to satisfy another identity 
that can be derived by considering the response of Tk on gauge transformations of the 
background field only. This background-field identity is in close relation to the modified 
Ward identity [pl| (for an explicit proof in QED, seej^TJ), and also imposes a constraint 
only on r| auge similar to Eq. fllOl) . As has been shown in fl9|], this identity does not cause 
further fine-tuning problems which would add to those that are posed by the modified 
Ward identity. 

In summary, solving the flow equation (P with the truncation (|5]) will result in an action 
functional rj^LA = A] which is invariant under the background-field transformation. By 
neglecting r| auge , this invariance is not identical to full quantum gauge invariance even at 
A = A, since the flow is not completely compatible with the modified Ward identities. 
This work is based on the assumption that these violations of quantum gauge invariance 
have little effect on the final result. 



3 Evaluation of the truncated flow 

We shall now solve the flow equation (|]) within the truncation ([5]) (neglecting r| auge ) and 
with T k nv as given in Eq. (|2|), 

TT[A]=!w k {6\ w k {e) = Y^e\ (11) 

where 6 := \F^ V F^ V . An important ingredient of the flow equation is the cutoff function 
Rk, which we display as 

x 

R k {x)=xr(y), y := — — , (12) 

Z k k z 

with r(y) being a dimensionless function of a dimensionless argument. We include wave- 
function renormalization constants Z k in the argument of r(y) for reasons to be discussed 
below. Note that Z k as well as R k itself are matrices in field space; different field variables 
may be accompanied by different Z k s and -R^'s. The cutoff function Rk has to satisfy the 
following standard constraints: 

lim Rk(x) > 0, lim Rk(x) = 0, lim Rk(x) — > oo, (13) 

x/k 2 ^0 k 2 /x~^0 k^A 

which guarantee that Rk provides for an infrared regularization, ensure that the regulator 
is removed in the limit k — > 0, and control the ultraviolet limit where = S\ should 

approach its initial condition S\ at the initial ultraviolet scale A. These constraints are met 
by the representation (0) and translate into constraints for r(y). Since we will identify 
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the argument x with the full T k at the background field, the first constraint of fll3|) must 
be formulated more strongly, 

lim R k (x)=Z k k 2 , r (y->0)->i (14) 

in order to guarantee that the one-loop approximation of the flow equation results in the 
true one- loop effective action. We shall not specify r(y) any further until we employ an 
exponential cutoff for the final quantitative computation (see. Eq. ( p.4|) ). 

Within the approximations mentioned above, the flow equation (|T5| ) can be written as 

d t T k [A = A,A] = -STr dtRk ^ 



2 if + i4(rf; 



isTr 



,(2) 



v)h(y) + ^-(g{y)-h{y)) 



k 



.(2) 



(15) 



where we abbreviated 



./ \ ~yr (y) ( , r(y) 

h (y) := T~s — FT' 9\y) ■= t- — rr- (16) 
l + r(y) l + r{y) 



In Eq. (pT5|), we also defined the anomalous dimension 



T] := -d t \nZ k = -^-d t Z k , (17) 

which is matrix valued in field space similarly to Z k ; different field variables can acquire 
different anomalous dimensions. We would like to draw attention to the appearance of the 
term ~ dtT k on the right-hand side of the flow equation. This term arises from writing 
T k into the argument of the cutoff function. It reflects the fact that the cutoff adjusts itself 
under the flow of the spectrum of T k 2 \^\ Now it is useful to introduce (at least formally) 
the Laplace transforms h(s) and g(s) of the functions h(y) and g(y): 

oo oo 

h(y) = [ dsh(s)e-*', g(y) = [ ds~g(s) e -y s . (18) 



These Laplace transforms h(s) and g(s) can be viewed as cutoff functions in Laplace space: 
they should drop off sufficiently fast for large s (small s) in order to regularize the infrared 
(ultraviolet). For instance, the infrared constraint (0) translates into 

oo oo 

h(0) = J dsh(s) = 1 = J dsg(s) = g(0). (19) 



2 Although Ifc was also used as the argument of the cutoff in [|l2|, the term ~ d t T k 2 ^ has been neglected 
in that calculation. The necessity of this term was pointed out to us by D.F. Litim. 
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Additional useful identitiesjor these functions are discussed in Appendix [D| Furthermore, 
introducing the functions H(s) and G(s) by 

±H(s) = h(s), H(0) = 0, 

±G(s)=g(s), 5(0) = 0, (20) 
a convenient form of the flow equation can be found which reads: 

d t T k = ±J^(H(s)-G(s))d t STrex V (-s^ k 



(2) 
k 





oo 



2 



+^ J ds~g(s)STr(2- V )expl-s^). (21 



The great advantage of this form is that the right-hand side of the flow equation has been 
transformed into a propertime representation.^ The (super-)trace calculation reduces to a 
computation of a heat-kernel trace for which there are powerful techniques available. 
Within the actual present truncation (neglecting rf auge ), 



T k [A, A) = I W k {9) + rf [A, A] + Tf[A, A], (22) 

J X 



(2) 

T k still has a complicated structure which inhibits obtaining general and exact results for 
the heat-kernel trace. Fortunately, a general solution is not necessary; we merely have to 
project the right-hand side onto the truncation, implying that we need only the dependence 
of the right-hand side on the invariant 9 = \F®F® \ other invariants occurring in the heat- 
kernel trace are of no importance in the present truncation. Now the crucial observation is 
that the heat-kernel dependence on 9 can be reconstructed by performing the computation 
for the special field configuration of a covariant constant magnetic field (as it is explicitly 
defined in Eq. ( |A.1| )). In addition to this, we can perform the computation for A = A. 
For this field configuration, the flow equation finally depends only on the field parameter 
^ = \F^ v F^ y = i-B 2 , where B denotes the strength of the magnetic field; the latter is 



3 This representation of the flow equation should not be confused with the so-called propertime RG p2| . 
The latter represents a RG flow equation that is derived by RG improving one-loop formulas in propertime 
representation, and has been used in a variety of studies p3| . However, a propertime flow is generally 
not exact, as was proved in p4j , j25|: generic propertime flows can neither be mapped onto exact flows 
in a derivative expansion nor correctly reproduce perturbation theory. By contrast, our flow equation is 
derived from an exact RG flow equation and corresponds to the generalized propertime flow proposed in 
p5| . The essential difference to (standard) propertime flows is the inclusion of dtT^ terms. In the present 
work, particularly these terms will be important and will not be neglected. In agreement with |25| , our 
findings therefore suggest that propertime flows may be improved towards exact RG flows by including 

(2) 

the ~ dt^k terms systematically. 



10 



pseudo-abelian and points into a single direction in color space, characterized by a color 
unit vector n a . Extracting this B dependence of the heat kernel allows us to reconstruct 
the flow of W k (\B 2 ). 

It should be stressed that considering a covariant constant background field is nothing 
but a technical trick to project onto the truncation; we do not at all assume that such a 
background represents the vacuum configuration of Yang-Mills theory, as is the case, e.g., 
in the Savvidy vacuum model |26| . 



(2) 

This trick also allows us to decompose the operator Y k into linearly independent 
pieces. In particular, the gauge-field fluctuations can be classified into modes with gener- 
alized transversal (T) and longitudinal (L) polarization with respect to the magnetic-field 
direction in spacetime, and into parallel || and perpendicular _L modes with respect to 
the field direction in color space. Introducing the corresponding projectors Pl,t and P\\,±_ 
(explicitly defined in App. (0)), the operator 1^ can be represented as 



,(2) 



-£> T 

a 



PtP±[W^V t ]+P l P ± 

+p t p,i [w' k {-d 2 ) + w k ' s] + p L P|| 



- (~d 2 ) 
a 



(23) 



where = [yl = A, A], and we drop the bars from now on. Here we also defined the 
operators 



(Pr)t = i-D%u + 2igF t 



\ab 



liv ) 



J fJ,U 



n a F%- 



(24) 



The formal symbol P g h in Eq. (|23|) projects onto the ghost sector, and W k = ^Wk(9); for 
details about this decomposition, see App. |A| 

At this point, we are free to choose different cutoff wave-function renormalizations 
for each of the linearly independent parts in Eq. (^). If we were solving the flow equation 
exactly, the final result would be independent of this choice; however, for a truncated flow, 
a clever choice can seriously improve the approximation. With regard to the form of Rk 
in Eq. ([12]) , it is obvious that the Z^s control the precise position at which the scale k 
cuts off the infrared of the momentum spectrum. Since the latter is determined by the 
Laplace-type operators — d 2 , — D 2 ,Vt in Eq. (p3|) , we can cut them off at k 2 by choosing 



J ghost, k 



'L,fc 



z 



a 



T,fc 



J F,k, 



(25) 



for ghost, longitudinal, and transversal fluctuations, respectively^ This choice guarantees 
that the longitudinal and ghost modes are cut off at the same point, providing for a 

In agreement with our approximation of neglecting r| augc , we also do not distinguish between the 
wave-function renormalization constant of the classical field A and that of the background field A; if r| augc 
were taken into account, this distinction would have to be made 27|. A similar problem would occur, if 
we relaxed the constraint (|l4|) for the cutoff; by keeping track of different wave-function renormalization 
constants, such a problem was explicitly solved in p8j. 
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necessary cancellation. As a side effect, the now becomes independent of the gauge- fixing 
parameter a, so that we can implicitly choose Landau gauge a — > 0, which is known to be a 
fixed point of the flow [[/]],|29]. Finally, the transversal cutoff wave- function renormalization 
is set equal to the gauge-field wave function renormalization, which can be read off at the 
weak-field limit, r fe L4]| w . f . ~ W£(0) = ^F^F^. 



Using trace identities found in flTj|, the heat-kernel trace occurring in Eq. can be 
further reduced to 



42) 



STr e 



Tr,,<' ^ 

+Tr a!cL e 



z F,k 



Tr xc e ^ " ' Tr,,.e "' . 



(26) 



where the traces can act on spacetime (x), color (c), or Lorentz (L) indices. For the trace 
in Eq. ( p3[) involving rj (matrix- valued) , all terms in Eq. (|26|) containing Zpj, will acquire 
an anomalous- dimension contribution which we will also call r\ for simplicity: 



-d t \nZ F<l 



1 



d t Z, 



(27) 



The various heat-kernel traces are computed in Appendix [B| In order to display the result 
concisely, let us define the auxiliary functions 



h{vx,v 2 ) 



1 



u d/2 
1 1 



(d-r 



sinh u 



2 u sinh u 



u 



2 u d l 2 sinh u ' 
1 . 

(1 — v 2 ). 



(28) 



Equipped with these abbreviations, the flow equation can be written as 

nv 2 -] 

/ rl$ < n( z\ 

2(4tt) 

o v L l — L 



d t W k {9) 



00 ( nv 2 -i 

n v L I — 1 



— (2 - 77) /; ' 



l[H(s)-G(s)]d t 



N 2 ~l 



d/2 
I 



2/3 



1=1 



d/2 



(29) 



where B\ = g\ui\B, g denotes the bare coupling, and v\ represents the I = 1, . . . , N 2 — 1 
eigenvalues of the color matrix (n a T a ) bc . The auxiliary functions f\ in the last line are 
understood to have the same arguments as in the first lines. It is convenient to express the 
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flow equation in terms of dimensionless renormalized quantities, 



9 



l.d-4 7 -\ -2 



= g 2 k- d z Fjk e = k- 4 g 2 e, 

M#) = 9 2 k~ d W k (6) = k- 4 Z~lg 2 W k (k 4 $/g 2 ) 



(30) 



and evaluate the derivative d t from now on at fixed d instead of fixed 9. As a result, the 



flow equation fl29|) turns into 



-(4 -rj) w k + Adw k {$) 



2(47r) d / 2 



ds h(s) 



N 2 -l 



^(/lW)-/2W)f-2/3 



3\SW k , ~ 



w k + 2$w k 



N 2 -l 



1=1 



-77 J dsg(s) 



oo 

+ / ds(h{s) -g{s)) 



2 ^ fi{sw k b{) bf 2 - f 3 ( sw k , - 



w k 



w k + 2i3w k 



(31) 



d t w k 4$w k 



w k 



w k 



N 2 -l 



2 fi( sw kbi) b'i 



d/2 



1=1 



-hd t - A^) f 3 ( sw k , - — ™ k „ 
d \ w k + ZifWk 



where Wkifi) = d^w k {'&)- l and we abbreviated bi = \vi\y2 , d. Equation (pT|) represents one 
of the main results of the present work. Within the chosen truncation, this flow equation 
leads to the full quantum effective action of Yang-Mills theory upon integration from its 
initial condition at A down to k = 0. 

As a first comment, we would like to mention that we rediscover the flow equation 
of P"2| if we perform an expansion for weak magnetic field and if we neglect all terms 



proportional to dfT k . In order to isolate the latter from the rest of Eq. (|3T|), the single 
factor of g(s) in the second line should be represented as h(s) — (h(s) — g(s j), and then all 
terms proportional to (h(s) — g(s)) should be dropped (cf. Eq. flT5|)). 

Obviously, the dtT^ terms ~ (h(s) — g(s)) modify the flow equation extensively.^ 
They seriously increase the degree of complexity of this partial differential equation, so 
that neither an analytic nor a numeric evaluation is straightforward. The next section will 
be devoted to a search for the simplest possible and consistent approximation. 

Finally, we remark that the flow equation contains a seeming divergence: in the limit of 
small k, the s integrand may not be bounded for s — > oo, owing to the last term ~ sinh u in 
the auxiliary function j\ given in Eq. (|28|). However, this divergence is well understood and 



5 Incidentally, it is easy to show that no admissible cutoff shape function r(y) exists such that h(y) 

(2) 

g(y). Hence, the dtT k terms are present for all cutoff shape functions. 
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can be controlled. It arises from the Nielsen- Olesen unstable mode [3(J in the operator X>t, 
and can be traced back to the fact that the gluon-spin coupling to the constant magnetic 
field can lower its energy below zero. Because of this mode, the covariant constant magnetic 
field is known to be unstable, if considered as the quantum vacuum state of Yang-Mills 
theory. The divergence can be identified as a pole at complex infinity. The s integral can be 
properly defined by analytic continuation, resulting in a real part as well as an imaginary 
part. The real part is indeed important because it contributes to the f3 function and the 
form of the effective action in a perturbative computation (see below). The imaginary part 
is interpreted as a measure for the instability of the constant-field vacuum. 

As we have stressed before, the constant-magnetic-field background is just a calcula- 
tional tool in the present context, and the validity of the flow equation is not based on this 
background. Therefore, the s integral can be properly defined by analytic continuation 
around this pole at complex infinity. The resulting real part will be a valid and important 
contribution to the flow, but the imaginary part is of no relevance here. If we were really 
interested in a constant-field vacuum, the flow generated by this imaginary part would de- 
scribe how the instability develops upon integrating out the unstable mode in a Wilsonian 
sense. 



4 Running gauge coupling in d = 4 

In order to find a strategy for solving the flow equation ([H]) within a first simple approxi- 
mation, let us take a closer look at the standard procedure employed for ordinary cutoffs 
without dtT k terms. In such a case, the partial differential equation can be rewritten 
as an infinite set of coupled ordinary first-order differential equations by expanding the 
truncation, e.g., 

oo 
i=l 

Note that, owing to the choice ( p5|) for and the definition ([H]), W\ — 1 is fixed, so 
that the generalized coupling W\ is traded for the anomalous dimension rj. As a result, we 
obtain infinitely many flow equations for the couplings Wi which, for an ordinary cutoff, 
read 

^ordinary cutoff = W 2, • • • , W i+1 ), i = 2, 3, . . . . (33) 

Equation (^) is supplemented by an additional equation for rj. The functions are 
obtained as the zth coefficient of the •& series expansion of the flow equation's right-hand 
side. This infinite tower of equations is then approximated by a finite one by setting all 
Wi = by hand for some i > itrunc, resulting in itrunc equations for itrunc variables. The 
quality of this further truncation can be checked by varying itrunc- 

This recipe cannot be directly applied to the present case involving the spectrally 
adjusted cutoff because an expansion of the flow equation ( |3~TD will be of the form 

d t Wi = Xi(rj, w 2 ,..., w i+ x) + Yi(rj, w 2 , ■ ■ ■ , w i+1 ; d t w 2 , • • • , d t w i+1 ), i = 2, 3, . . . . (34) 
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It is tempting to truncate this tower by setting not only u>i>i trunc = by hand, but also 
(9tW7i>j trunc = 0. This is too naive, however, because all d t Wi, if understood as the left-hand 
side of Eq. (|34"D, receive nonzero contributions on the right-hand side, even if i > Ztrunc- 
Neglecting these right-hand sides would correspond to neglecting some w^s which are in 
the truncation i < Ztrunc- 

In order to apply the above-mentioned recipe, we have to bring Eq. (|34|) into the form 
of Eq. (|33|), i.e., we have to solve for the d t w^s. Formally, this is possible by observing 
that the functions Yj, as they are derived from Eq. (pT|), are linear in all d t Wi and rj, and 
the Xi are also linear in rj. Introducing a "vector" w t with components 

^ := { Zl = d t l t for i = 2,3,... }' (35) 

equation fl3H]) can be written as0 

w u = Xi(w 2 , ■ • ■ , Wi+i) + Y ij (w 2 , • • • , w i+1 ) w tj , (36) 

or symbolically, w t = X + Y -w t . Provided that the operator 1 — Y is invertible, the desired 
solution is formally given by 

wt = -^y ■ X, (37) 

where the right-hand side is a function of w 2 , W3, . . . only. Now, the approximation strategy 
for the ordinary cutoff can be applied to Eq. fl3~71). Nevertheless, the resulting finite tower of 
differential equations is substantially different from the ordinary case, even for the smallest 
^trunc- This is because Xi and Yij are generally nonzero, even for i,j > i t runc, since they 
depend on the remaining Wi<i tiunc (and numbers such as d and N). And since they are 
infinite dimensional, we find an infinite number of terms on the right-hand side of the flow 
equations, in contrast to a finite number for ordinary cutoffs. 

For the remainder of this section, we shall evaluate Eq. (|37|) in the simplest possible 
way by neglecting all w^s with i — 2, 3, . . . and retaining only the anomalous dimension 
rj, which is related to the (3 function of Yang-Mills theory via 

P(g 2 ) = d t g 2 =(d-4 + r])g 2 , (38) 

so that in d = 4 we simply have f3(g 2 ) = f]g 2 . We would like to stress that the approxima- 
tion of neglecting all w^s at this stage is not at all equal to neglecting them right from the 
beginning. This further truncation is only consistent after we have disentangled the flows 
of all WiS by virtue of Eq. ( |37|) . 

For the investigation of the rj equation, corresponding to the first component of the 
vector equation fl37D , it is useful to scale out the coupling constant, so that X and Y no 



6 The meaning of the quantities Xi and changes here slightly, because the 77 and dtWi dependence is 
pulled out compared to Eq. (|34|) . 
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longer depend on the coupling: 

X^GX, Y^GY, G:=^£yj- 2 . (39) 
In d = 4, the convenient coupling G is related to the standard strong coupling constant 

2 

a s = |- = 87rG. Using Eq. (P§|), we can perform a "perturbative" expansion of the r\ 
equation: 

-7/ = u> fl = ( y^gy ) GX 3 = G(l + GY + G 2 Y 2 + ...) l3 X ] 

= G ( G rn Y m j Jfy (40) 

\m=0 / y 

The explicit representation of Y and X can be found by inserting the expansions developed 
in Appendix |C| into Eq. fl3T|), and performing the propertime s integration; the latter results 
in the moments hj,gj of the cutoff functions h(s),g(s), 

oo oo 

hj := J dss^h^s), g 3 := J dss j g(s), (41) 
o o 

which are discussed in Appendix |B|. In conclusion, we find: 



X % = -¥ +1 r t h 2t _ d/2 i\ [ (d-2) y "." J B 2i 



(2 2i - 2) 



(2i)! (2z-l)!. 
Yij = Aij + -B^ + Cij, (42) 

where i? 2 i denotes the Bernoulli numbers, and the auxiliary matrices A,B,C are given by 



i,j = 1,2, 


...): 


... 


A^ 

1 











2 n r n (/l 2 n-d/2 - 92n-d/2) ({d - l) 2 j^B 2n - (2^TT)! 



- n=l+i—j 



V (2i)l (2i-l) 

9* = 0i£ /.t i + 1 , ■ (43) 



These explicit representations ( |42"D and (^) can be inserted into Eq. ([40]) , and the anoma- 
lous dimension and the /? function can be computed straightforwardly to any finite order 
in perturbation theory within our truncation. 
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As an example, let us compute the two-loop (3 function in d = 4 spacetime dimensions 
for SU(iV) gauge theories: 



2 2 22iV <? 4 /77iV 2 127(3iV 2 -2) ^ , \ <? 6 

= % =-^(4^"^^ 45 (^-^)^J(4^ + - • 

(44) 

Here we already used the fact that ho = 1 = go are independent of the shape of the 
cutoff function, so that the one-loop coefficient turns out to be universal as it should be 
and agrees with the standard perturbative result; this should serve as a (rather trivial) 
check of our computation. Within our truncation, the two-loop coefficient does depend 
on the cutoff function. In order to compare with [12j where the dtT^ terms have been 



neglected, we choose the exponential cutoff defined in Eq. ( p.4| ), implying that g_ 2 = 1, 
h_ 2 — 2<C(3) — 2.404, and h 2 = 1/6 as computed in Appendix |D]. From Appendix [E|, we 
take over that t^ =2 = 2 and t% =z = 9/4. Inserting all these numbers and comparing this 
to the perturbative two- loop result, 

a . 2 , 22N g 4 68N 2 g 6 

/W(S) = -— 7^- — 7^ + -.., (45) 

we find a remarkable agreement of 99% for the two-loop coefficient for SU(2), and 95% 
for SU(3). This should be compared to only a 113% agreement of these coefficients in the 
case where the d t V k terms are neglected p2| . The inclusion of these terms appears to 



represent a serious improvement. 

However, the picture is not so rosy as it seems to be in view of this result. The reason 
is that our two- loop result is cutoff-scheme dependent, and we may easily choose a cutoff 
with a worse agreement at two loop.Q Only recently has it been explicitly shown how 
to obtain the correct scheme-independent two-loop (3 function within the framework of 
the exact RG ||27|; for this, a careful distinction has to be drawn between the running of 



the coupling with respect to k or the RG scale \x (see also pip. Within our truncation 
here, we can nevertheless turn the argument around by remarking that the exponential 
cutoff is obviously well suited for the present truncation in the sense that it minimizes the 
combined effect of the neglected terms such as u> 2 , (F° U F£ V ) 2 , or r| auge , etc. on the two-loop 
coefficient. This particularly justifies the use of the exponential cutoff for an investigation 
of the complete sum in Eq. flip]) and the strong-coupling domain. 

Let us summarize what has been achieved so far: in order to extract the flow of the gauge 
coupling, the flow equation ([31] ) has to be studied near $ = 0, where the information about 
7] is encoded. This suggests an expansion of Wk{"d) in powers of leading to completely 



disentangled flow equations (|37 ) for the generalized couplings 77, w%, W3, .... However, since 



7 Actually, this point is even more subtle |27j: cutoff-scheme independence of the two-loop /3 function 
coefficient holds only for mass-independent regulators. The regulator obviously does not belong to this 
class, so that cutoff-scheme dependence has to be expected. However, since we have no information about 
the true two-loop coefficient for our /c-dependent regulator Rk, we shall still use the two-loop coefficient 
for a mass- independent scheme as our benchmark test. 
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the original flow equation fl3T| ) is represented as a parameter integral, its expansion can be 
asymptotic, which implies that the series expansion in terms of the coupling G in Eq. ( flUD 
will be asymptotic as well. This agrees with the general expectation that perturbative 
expansions of quantum field theories are generically asymptotic. 

In practice, this means that the coefficients (for later convenience, we shift the index 
m here) 

a m := -(r- 1 )^, m = l,2,... (46) 

in Eq. ( flOf ) grow rapidly, so that any arbitrarily large but finite truncation of the series does 
not make sense. It turns out that these coefficients grow even more strongly than factorially 
and alternate in sign for the exponential cutoff (we shall comment on other cutoffs later). 
This does not mean that any physical meaning is lost, but, loosely speaking, that we have 
expanded an integrand which we should not have expanded. Yet, there are well-defined 
mathematical tools for reconstructing the integrand representation out of the diverging 
sum [32]. In other words, we are looking for a (well-defined) integral representation that 
upon asymptotic expansion leads to a series that agrees with Eq. fl4"0p. As is known also 
from various physical examples |33[], just taking only the leading growth of the coefficients 



into account leads to a good approximation of the integral representation. 

Concerning the coefficients a m , the leading growth (l.g.) can be isolated in the term 
that contains the highest component of X, i.e., X m , yielding 

<h s ' = —Xi, a 2 ' g ' = —Y 12 X 2 , a^f = — ^12^23^34 • • • Y m _ l m X m . (47) 
Inserting the representations fl4"3|) into Eq. ( |4"?D for the exponential cutoff, we find 

.1* „ „_>,„-ir(™ + 3(iv 2 -i))r( m . + i)_ f n 2 2m -2„ 4 



* = ^r~^ '- (2 - pjssj ) ■ <«) 

where we abbreviated c = 2£(3) — I. Let us first concentrate on SU(2), where r m = 2 for 
m = 1,2, . . . (see Appendix |Ej); let us nevertheless retain the N dependence in all other 
terms in order to facilitate the generalization to SU(3). 

Actually, Eq. (flST) also contains subleading terms. First, we observe that the last term 
~ l/r(2m) is negligible compared to the term ~ B 2m for large m. Nevertheless we also 
retain this subleading term, since the m = 1 term contributes significantly to the one-loop 
/3-function coefficient which we want to maintain in our approximation. Furthermore using 
the identity 

(-l)"T(2m-l) . 
B ^2 = 2 L- 3 \ 2m -2 C(2m - 2), (49) 

it is tempting to use the ^-function representation 



00 1 



<-i l 
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and retain only the / = 1 term, since the others are subleading. Whereas this approximation 
is indeed justified for the Bernoulli number B 2m in Eq. (^8j), we have to retain the full ( 
function for the -B2771-2 factor, since here we encounter the ( function at zero argument for 
m = 1 where Eq. (^) is no longer valid. In conclusion, we resum the complete coefficient 
displayed in Eq. ([IB]), including the leading and also subleading terms. 
In the spirit of Borel summation f3~2fl , we introduce integral representations for the 
special functions occurring in Eq. (|48"1), in particular the representation [j34H 



r(2m - 1) C(2m - 2) 



1 



1-2 3 - 



-2m 



dt- 



e* + 1 



t 



2m-2 



m > 1/2 



(51) 



for the ( function in Eq. ([4£]), and an Euler B function representation for a combination 
involving the last term in Eq. ([48]): 



B(m, m) 



r(m + 3(iV 2 - l))r(m+ 1) _ T(m + 3(N 2 -1)) T(m + 1) 



T(2m) 



T(m) 



T(m) 



m(m + 1) . . . (m + 3iV 2 - 4) • m ds s m ~ L (1 - s 



\m— 1 



ds 



— ^ 

ds J 



(37V 2 -3) 



ds' 



(52) 



s'=l-s 



For the remaining T functions, we use the standard Euler representation. Exploiting these 
identities, we are able to resum Eq. ([40]) to this order: 



I.e. nin 



g ( (: '" 1 y '" 1 ) x j - Yl a -' 

\j m=l 



\m=l 
V* + Vb, 



(53) 



where r] a is related to the term ~ 5 2m in Eq. (f4*B[), whereas t]^ is related to the term 
~ l/r(2m). The integral representation of r/ a reads 



32NG 



e*-(»i+- a ) SlS 3JV 2 -3 



5 



cGs\S2t^ 

2tt 4 / 2 



2^ 



cGs\S2t* 

472 



(54) 



where we defined the sum 



~ .Z^ 1 _ 23-2m 



m=l 



2i + X 

i=0 + 2* 



=: 5(g). 



(55) 
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The first sum arises from the asymptotic expansion and is strictly valid only for \q\ < 1; 
however, the second sum is valid for arbitrary q, apart from simple poles at q = — 2 2j , and 
rapidly converging, so that this equation should be read from right to left. 

The second part 7]^ deserves a comment: as it arises from the last term in Eq. (|48|), 
~ l/r(2m), it originates in the last term sinhw of the auxiliary function fx in Eq. (|28|), 
which stems from the lower end of the spectrum; in particular, it contains the Nielsen- 
Olesen unstable mode. This mode is reflected in a simple pole in the following integral 
representation for rft>- This pole gives rise to an imaginary part of the full integrand. As we 
have stressed above, the imaginary part created by this unstable mode is of no relevance for 
the flow equation here, so that the proper treatment of the integral results in a principal- 
value prescription maintaining the important real part. For a numerical realization, this 
prescription can best be established by rotating the t integral arising from Eq. (|5T| ) by an 
angle of, e.g., 7r/4 from the real axis into the upper complex plane and then taking the real 
part. In conclusion, we get: 



N _ 2 32NG „ T, , e^* fdf N ~ 3) d w2 , , n .cGss't 2 

Vh = ; 5 rRe Ids l^-dt^rn — — s 3N ~ 3 s'S - — 

/b T 3iV 2 -2 J J * ( e ^*+l) 2 W ds ' V 27r2 



(56) 



Although it seems that we have seriously complicated the problem by trading the single m 
sum in Eq. ( PH) for a number of integrals and sums, we stress that all integrals and sums 
in Eqs. (|54|) and (]5"6|) are finite and well defined. 

Before we present numerical evaluations of these integrals and sums, let us discuss 
some features analytically. For small coupling G = g 2 /[2(47r) 2 ], we can again expand the 
integrals asymptotically and obtain 



»° = l N (h? + ---' ,h = - m jtf- + -' (57) 

so that we rediscover the one- loop (3 function (cf. Eq. (0)) as a check. Next, we observe 
that ?7 a (containing the true leading-order growth of the o^'s) is positive not only for small 
but arbitrary G. In order to extract large-G information, we note that the sum S(q) can 
be fitted by 



S(q)~c lx l^-c 2 ^ + q, Cl ~ 2.27, c 2 ~ 0.7 (58) 

within 1% accuracy, implying that S(q) ~ C\ y/q for large q. In this limit, which corresponds 
to large G, r] a can be evaluated analytically and we find 



/f C ( 3 ) ln2r(5/2) r(3iV ' 3/2) G 3/2 ~ 3.24G 3 / 2 . 
V 2 SV ; v 1 ' T(3N 2 - 2) 
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Figure 2: Anomalous dimension rj = /3(g 2 )/g 2 for SU(2) Yang-Mills theory in d — 4 versus 
G = g 2 /[2(47r) 2 ]. The long-dashed line represents the contribution r] a , the short-dashed 
line rj b , as defined in Eqs. fl53p and (|5ED; the solid line is the sum of both. 



Without going into details, we note that there exist fits for S(— iq) similarly to Eq. ( |5SD 
involving a square-root behavior for large q (large G). It turns out that the G 3//2 coefficient 
vanishes exactly, so thatQ 

^=2 „ +G l/2. (60) 

Obviously, r/b is subleading for large G, which agrees with the fact that it arises from 
subleading parts in the coefficients a^f . Moreover, rjb becomes positive for large G, so that 
there should be a zero in between. This is already the first sign of an infrared stable fixed 
point at which r](G*) = 0. 

For a numerical evaluation of r] a and r/b, we employ the representations given in Ap- 
pendix |F|. We depict the anomalous dimension rj and its parts r] a and rjb in Fig. [| for the 
gauge group SU(2). The plots agree with the analytical estimates given above, and we find 
an infrared stable fixed point at 

Gf =2 ~0.45 =► «f= 2 ~ 11.3. (61) 

By virtue of Eq. (|38|), the running gauge coupling approaches this fixed point upon lowering 
the scale k in the infrared, implying scale invariance. The complete flow of the coupling is 
obtained by integrating (3{g 2 ) = d t g 2 = r\g 2 and has been plotted already in Fig. 0. 

For the gauge group SU(3), we do not have the explicit representation of the color 
factors r m at our disposal. As discussed in Appendix [FJ, we instead study the two extremal 
cases for the color vector n a pointing into the 3 or 8 direction in color space. Inserting 

8 We do not evaluate the precise coefficient of the G 1//2 here, since the large-G expansion introduces 
artificial singularities for the t integration at t — > 0. A more careful treatment reveals that the coefficient 
is positive. 
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Figure 3: Anomalous dimension rj = /3(g 2 )/g 2 for SU(3) Yang-Mills theory in d — 4 versus 
G = g 2 /[2(Att) 2 ]. The black lines correspond to rj^ =3 , the grey lines to ti^ =3 , as defined in 
Eq. (|S2"D. The meaning of the dashed lines is as in Fig. |^. 



the corresponding quantities r/3 3 or r/g 3 as found in Eq. ( E.5|) into Eq. (|^) allows us to 
display the anomalous dimension r/ Ar=3 in terms of the formulas deduced for SU(2): 



N=3 ^ N=2 



+ V= 2 

AT-+3 3 



N=3 N=2 
V 8 = V 



N^3,c^3c/4 



N->3,c->c/4 

(62) 



The notation here indicates that the quantities N and c = 2((3) — 1 appearing on the 
right-hand sides of Eqs. (^4j) and (|56"D will be replaced in the prescribed way. Figure || 
depicts our numerical results, and we identify the position of the infrared fixed point in 
the interval 

G? = 3 = [^,8, G %3 ] ~ [0.225, 0.385] => a? = 3 ~ [5.7, 9.7]. (63) 

This uncertainty of the precise position of the fixed point is not a shortcoming of the 
techniques involved (e.g., using the covariant-constant magnetic background), but is due 
to our ignorance of the exact color factors r m . 

Let us conclude this section with some remarks on the resummation: first, we should 
stress that the results for the fixed point are derived from a resummation of leading and 
subleading parts of the complete asymptotic series (pEOf ). We have checked that the sub- 
subleading parts (not included in the present resummation) alternate in sign, so that their 
contribution will be regular. However, we were not able to systematize the sub-subleading 
terms in a way that they can be resummed further in a consistent way. 

Secondly, we performed the computation for the exponential cutoff. It would be de- 
sirable to test the stability of the fixed point by using different cutoffs. Unfortunately, 
we could not find another cutoff shape function r(y) for which the resummation could be 
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done. For many cutoff shape functions used in the literature, the series is also asymptotic 
and alternating, but the sign changes not from one coefficient to the other but from one 
group of coefficients to the next, i.e., a ni . . . a n2 > and a n2+ i . . . a n3 < for n\ < n 2 < n 3 , 
etc. The outstanding role of the exponential cutoff may be attributed to its close relation 
to the Bernoulli numbers and their properties. 

Let us finally stress once more that the generalized Borel resummation of Eq. fl40|) does 
not represent an uncontrolled extrapolation of finite-order perturbation theory. As we have 
the exact all-order result at our disposal, the resummation corresponds simply to an - also 
mathematically - well-defined transformation of a series into an integral. 



5 The role of the spectrally adjusted cutoff 

This short section is devoted to a heuristic discussion of the special role played by the 
spectrally adjusted cutoff in this work, focusing on the truncation employed. 

The spectral adjustment of the cutoff function to the spectral flow of T k arises from two 

(2) 

sources: first, from using T k in the argument of Rk and, second, from including a carefully 
chosen wave-function renormalization constant in the cutoff. The latter technique is 
well known in the problem of calculating anomalous dimensions in scalar and fermionic 
theories. In order to get a feeling for these two improvements, let us first consider the flow 
equation, neglecting all d t T k 2 ^ and d t Zk terms. For the anomalous dimension, we would 
then obtain: 

v = - bo ji^- bi jk~ 2W2 ' (64) 

with some coefficient b\, and bo being the correct one-loop result. Obviously, choosing 
furthermore the truncation w 2 ,Ws, ■■■ = leaves us with a purely perturbative lowest- 
order result. This means that all nonperturbative information is contained in the flow of 
w 2 , which in turn can be reliably computed only by including w 3 , etc. A good estimate 
therefore probably requires a very large truncation. Even if the precise infrared values of 
the higher couplings W{ may not be very important, their flow exerts a strong influence on 
the running coupling in this approximation. 

Let us now take the d t Z^ terms into account, but still neglect the d t T k 2 ^ terms. In this 
case, the flow equation results in the following expression for the anomalous dimension: 

„_ 6 W + 61 (£r w * (65) 



1 + di (^2 + d 2 J £ f w 2 



with further coefficients di,d 2 , where d± < 0. Particularly this d\ makes an important 
contribution to the two-loop (3 function coefficient. Contrary to Eq. (0), this equation 
contains information to all orders in g 2 , even for the strict truncation w 2 ,w 3 , • ■ ■ = 0. We 
have to conclude that an adjustment of the cutoff function using a cutoff wave-function 
renormalization Zk is an effective way to put essential information of the flow of the higher 
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couplings W2,w?,. . . into r/. In other words, the truncated RG trajectory better exploits 
the degrees of freedom left in the truncation. Let us note in passing that the flow governed 
by Eq. (|65|) runs into a kind of Landau pole for ~ 1 / 1 c?i | , even if the flows of w 2 and 
higher couplings are included. This "disease" has occurred in many flow equation studies 
in Yang-Mills theory [|, §, Q, §5| . 



(2) 

Now let us turn to the full flow equation, including the terms generated by d t T k . As 
explained in the previous section, the right-hand side cannot be displayed in terms of the 
wfs in closed form, because infinitely many terms contribute. Even if we set all wfs to 
zero, which indeed corresponds to our final approximation, the anomalous dimension reads 

_ b 0J £ 7 + b 1J £^ + b 2J £ 1 , + ... 

1 + di + d 2 (4^4 + d 3 ^ + . . . 

with some real coefficients bi and d^ (expanding Eq. (J^) in powers of g 2 results in Eq. ([40D ). 
Whereas the nonperturbative dependence of 77 in Eq. ( j65| ) resembles that of a Dyson series 
and is controlled by one coefficient (d\ in that case), Eq. ( |66|) contains nonperturbative 
information from infinitely many coefficients. The latter arises from the flows d t Wi which 
all contribute to Eq. (|66|). We conclude that the spectrally adjusted cutoff provides for 
an efficient reorganization of the flow equation, so that a small truncation can contain 
information which, for ordinary cutoffs, is distributed over infinitely many couplings of 
a larger truncation. From this observation, we conjecture that the spectrally adjusted 
cutoff selects a truncated RG trajectory which is "optimized" with respect to the degrees 
of freedom within a chosen truncation. We furthermore conjecture that this trajectory 
does not flow into regions of theory space where the exact flow would be mainly driven by 
couplings which are not contained in the truncation, but is always driven by the couplings 
within the truncation in an optimized way. Whether the truncated RG trajectory flows to 
the true quantum action or not depends, of course, on the quality of the truncation. 

We have obviously verified these conjectures only for a truncation (if) within a trunca- 
tion (r],Wi). In fact, in order to exploit the properties of the spectrally adjusted cutoff, we 
first have to discuss the flow of a larger truncation, then disentangle the flows of the single 
couplings and finally restrict the calculation to the most relevant part under consideration. 

Let us finally point out that using the spectrally adjusted cutoff necessarily requires 

(2) 

introducing a background field, because T k in the cutoff function is not allowed to depend 
on the actual field variable. A background field generally complicates the formulation, and 
the technical advantages of the spectrally adjusted cutoff may be compensated for by these 
further complications. Gauge theories, however, may serve as a natural testing ground for 
the spectrally adjusted cutoff, since the background-field formalism is advantageous here 
for further reasons. 



6 Conclusions 

Starting from the exact renormalization group flow equation for the effective average ac- 
tion in SU(iV) Yang-Mills theory in d dimensions, we derive within a series of systematic 
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approximations the (3 function of the gauge coupling. In d = 4 spacetime dimensions, 
the resulting flow of the gauge coupling exhibits accurate perturbative behavior and ap- 
proaches a fixed point in the infrared. The fixed-point results are displayed in Eqs. fl6T|) 
and flUD for gauge groups SU(2) and SU(3). 

In view of the approximations involved, a number of improvements are desirable in order 
to confirm the existence of the infrared fixed point. Above all other possible improvements, 
such as enlarged truncations and explicit cutoff-shape independence (or insensitivity), a 
better control of gauge invariance under the flow is necessary. 

Nevertheless, in view of the flow equation studies performed in the literature for gauge 
theories so far, it is already remarkable that our approximation to the exact flow equation 
is integrable down to k — > 0; in many instances, the truncation revealed an explicit insuf- 
ficiency by developing a Landau-pole type of singularity at some finite k in one or more 
couplings. The new technique in the present work is the use of a cutoff function that ad- 
justs itself permanently to the actual spectrum under the flow. From a practical viewpoint, 
this cutoff condenses information, which is usually distributed over the flow equations of 
infinitely many couplings, into the flow equation of a single coupling (in this case the gauge 
coupling). We have reason to believe that the information, which is reorganized in this way 
into a single flow equation, is the relevant information that mainly drives the flow of the 
corresponding coupling. The fact that we improved the agreement with the perturbative 
two-loop running from merely 113% to 99% for SU(2) (using the exponential cutoff shape 
function) may serve as a hint in this direction. 

If the fixed point exists and our truncation even covers the true mechanism, it is still 
unlikely that our present results for a* are also quantitatively correct. We expect a low- 
ering of a* for larger truncations owing to the following argument: in our calculation, the 
position of a* is strongly governed by those modes which are also responsible for asymp- 
totic freedom (contained in rfo). If, in a larger truncation, operators of higher order are 
generated under the flow, these modes will generically lose influence, and the effects of 
the remaining spectrum contributing to r) a will be enhanced. This will shift a* to smaller 
values. 

A similar effect occurs upon the inclusion of quark degrees of freedom. The perturbative 
quark contribution to the f3 function is already positive. And since no ultraviolet stable 
fixed point is known in QED, we also do not expect negative quark contributions beyond 
the perturbative regime. Therefore, we expect not only the presence of the fixed point in 
full QCD, but also a substantial shift towards lower values of a*. Work in this direction is 
in progress. 

A comparison of our result with the literature is in order now, although it is generally 
difficult, owing to the various nonperturbative definitions of the gauge coupling; different 
definitions may agree perturbatively, but differ beyond perturbation theory. Our definition 
is standard in pure continuum gauge theory; moreover, it is equal to the interaction strength 
of static quarks with the gauge field. Nevertheless, it is not immediately clear to us how it 
can be related to a definition which is used, for instance, in lattice gauge theory p5 |. This 
may serve as a word of caution. 

The notion of an infrared fixed point for the gauge coupling has been used extensively 
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in recent years, especially in connection with the phenomenology of power corrections in 
QCD [|37|| . Furthermore, such a so-called freezing of the coupling has been discussed in 
phenomenological low-energy models [3S], and deduced from an analysis of the famous 
R e + e - ratio [j39fl . 

There are also various theoretical arguments favoring an infrared fixed point, e.g., 
even within a perturbative framework for a finite number of flavors Furthermore, 
investigating analyticity properties in the time-like and space-like (Euclidean) region, a 
scheme called analytic perturbation theory has been proposed, yielding an infrared finite 
coupling fill ; this program has been successfully applied to hadron and lepton-hadron 
phenomenology Having the above-mentioned reservations in mind concerning the 



various different nonperturbative definitions of the coupling, the question of how they are 
related to each other deserves further study. 

Moreover, an actual nonperturbative computation of gluon and ghost propagators has 
been set up in the framework of truncated Schwinger- Dyson equations in Landau gauge 
43j , revealing an infrared fixed point; these results also receive some support from lattice 



calculations MM. Again, the relation to our results is not immediately obvious, since the 



running coupling as defined in [|3| is obtained from the ghost-gluon vertex; furthermore, 
a nonperturbative treatment of the ghost sector turned out to be crucial in that work, 
but the four-gluon vertex was neglected. Nevertheless, there are also similarities: on very 
general grounds, it was found in the approximation of [||[] that the fixed point scales with 
the number of colors as a* ~ 1/N. We observe that the central value of our SU(3) result 
and the SU(2) result fulfil exactly this relation, although this is far from self-evident in our 
calculation. 

Let us finally discuss further implications of our result: comparing the full (3 function 
with its perturbative counterpart, we observe a quantitative agreement up to a s ~ 1. This 
does not, of course, justify the use of perturbation theory up to a s ~ 1 in general, but may 
explain why perturbation theory gives an accurate answer to some questions, even at its 
validity limit. 

Concerning the low-energy fixed-point region, one may ask whether our result provides 
for some signals of confinement and an expected mass gap in gauge theories. In the first 
place, the answer is no, since a strong coupling does not necessarily imply confinement. 
It is rather likely that the strong coupling of the gauge fields is necessary to give rise to 
a change of the effective degrees of freedom. These degrees of freedom (not necessarily 
included in our truncation) with probably nontrivial topological properties will then act 
as "confiners" . Also the picture of confinement arising in the framework of Landau-gauge 
Dyson-Schwinger equations ||3| cannot be contained in our truncation, since it is based on 
an infrared enhancement of the ghosts which are treated rather poorly in the present work. 
Improvements in this direction are also subject to future work. As far as a mass gap is 
concerned, the infrared fixed point behavior is compatible with such a gap; this is because 
a mass gap cuts off all quantum fluctuations of lower momentum, so that nothing remains 
to drive the flow. But the mere existence of an infrared fixed point does not require a mass 
gap. 

An indirect signal of a mass gap may be found in the analysis of the different spectral 
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contributions; as we have mentioned above, the perturbative (3 function is mainly deter- 
mined by the lowest modes in the spectrum, i.e., the lowest Landau levels in the covariant- 
constant field analysis. As is familiar from QED calculations, the lowest-Landau-level 
approximation is always appropriate if the field strength exceeds the mass of the fluctu- 
ating particle. This is certainly the case in the perturbative domain where the gluon is 
massless; hence the picture is complete. When we enter the infrared fixed-point region, 
the contributions from the remaining part of the spectrum r) a become important. In the 
Landau-level picture, this is always the case if a mass of the order of the lowest Landau 
level and beyond is present. The value of the mass then controls the influence of the re- 
maining spectrum. Therefore, the influence of the complete spectrum at the fixed point 
may be a hint for a hidden new mass scale in low-energy Yang-Mills theory. 



Appendices 

A Decomposition of 



Here we briefly describe the method developed in [12| for decomposing T k into smaller 
building blocks suitable for further diagonalization. The method is based on the observation 
that it is sufficient to consider only a covariant constant magnetic background field in order 
to project the flow equation onto the present truncation. 

The method consists of identifying those components of the quantum fluctuations which 
are appropriately oriented with respect to the background field; the latter is chosen to be 
of the type 

Al = n a A M , F M „ = d„A v - d u A^ = Be^ = const., (A.l) 

where n a is a constant unit vector in color space, n 2 = 1, and A M , F ^ denote the "abelian" 
gauge field and field strength. The constant tensor eh v characterizes the space directions 
which are affected by the constant magnetic field upon the Lorentz force, e.g., e^ 2 = — = 
1 for B pointing into the 3 direction. 

Let us first define two important operators involving the covariant derivative (D tJi [A\) ab = 
d tl 5 ab — igA c fl (T c ) ab in the adjoint representation: 

(Pr)% = (-D\ u + 2igF, u r b 

(V h )% = -{D®D)% = -D™D<?, (A.2) 



where (F^ iu ) ab = F^ v {T c ) ab . For covariant constant fields of the type (?Q) satisfying the 



equations of motion [D^, F^ v ] = 0, the operators and V L commute. As a consequence, 
projection operators can be introduced: 

P L = V^V L , P T = 1 - Pl, (A.3) 
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which obey P| 



TL — Pt,l, Pt + Pl = 1? -Pt-Pl = = Pl-Pt- The subscripts indicate that 
these projectors reduce to the standard longitudinal and transverse projectors in the limit 

Another pair of projectors can be defined which act solely in color space: 

Pf = n a n\ Pf = 5 ab - n a n b . (A.4) 

These four projectors are remarkably efficient in the present case; differentiating our trun- 
cation for Tk[A, A], as given in Eq. (|22|) , twice with respect to A and the ghost fields, then 
setting A = A and dropping the bar, we can represent the result as 



if [A 4 



1 



a 



PtP± [W^V t ]+P l P ± 

+P T P,I [W' k (-d 2 ) + W'i S] + P L P|, 
+Pgh [~D% 



a 



where we introduced 



(A.5) 



(A.6) 



and P g h projects trivially onto the ghost sector. 

Equation ( |A.5|) is perfectly suited for further manipulation, since the spectra of the 
operators occurring in the square brackets is known. This decomposition also offers the 
possibility of conveniently implementing different wave-function renormalization constants 
for each subcomponent. 



B Heat-kernel computations 

In this appendix, we summarize the results for the heat-kernel traces appearing in Eq. (|26|). 
Again, it is sufficient to perform the calculation for a covariant constant background field 
in order to disentangle the contributions to the flow of different operators. 

Let us first mention that all color traces occurring in Eq. (|26|) are of the form 

2V 2 -1 

tr c f(n c (TT b ) = Yl /M' (B-l) 
i=i 

where / is an arbitrary function, and v\ denotes the eigenvalues of the matrix (n c T c ) ab . 

We begin with the heat-kernel trace involving the Laplacian in the covariant constant 
magnetic background; the spectrum is given by 

Spect.P 2 : q 2 + (2n + l)B h Bi = g\vi\B, n = 0, 1, . . . , (B.2) 
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where g M denotes the (d — 2) dimensional Fourier momentum in those spacetime directions 
which are not affected by the magnetic field. The index n labels the Landau levels; their 
corresponding density of states is Bi/(2ir). Tracing over the spectrum, we obtain 

-Tr e-^-^-V* - - XBl (B3) 

n x ° ' 2(4vr)<V 2 \ d / 2 sinhA^' 1 ' 



Here, f2 denotes the spacetime volume. With reference to Eq. (fffi), the parameter A can 
be identified with A = sW k /(Z k k 2 ) or A = s/k 2 . 

Next, we turn to the heat-kernel trace involving the operator T>t as defined in Eq. ( A. 2 ). 
The spectrum is given by 

q 2 + (2n + l)Bi, multiplicity (d - 2) 

Spect.P T : q 2 + (2n + 3)B h multiplicity 1 (B.4) 
q 2 + (2n — l)Bi, multiplicity 1, 

with q and n as in Eq. (|B.2|) . The last line contains the Nielsen- Olesen unstable mode for 
n = |30 |, which has a tachyonic part for small momenta q 2 . Tracing over the spectrum, 
we find 

i l^xcL = E ^ ^ (d + sinh XB.) . (B.5) 

Finally, we need the following traces 

-Tr x e~^ - — - 



Q x 2(4vr) rf / 2 A rf / 2 ' 

iTr L e->^ 2 )-Vs = 2(d-l) J_ _2 1 X__ 

Q xcL 2(4vr) d / 2 \ d / 2 + 2(4vr) d / 2 \ d / 2 A + B 2 X' ' 1 J 



where S has been defined in Eq. ( |A.6| ) . Here and in Eq. ( |B.5| ) , the A parameters abbreviate 
A = sW' k /(Z k k 2 ) and A' = sW^/\z k k 2 ). 

Equations ( |B.3| ), ( |B.5| ), ( |B.6| ) serve as the main input for evaluating the right-hand side 
of the flow equation in Sect. [3]. 



C Expansions 

Here we shall explicitly display the expansions which are required for the analysis of the 
anomalous dimension in Sect. ^. The series given below are expanded in terms of the 
renormalized dimensionless field strength squared but they are also related to expansions 
in terms of the propertime parameter s or the renormalized coupling g 2 . Since we are 
expanding an integrand and then interchange integration with expansion, the resulting 
series can (and will) be asymptotic, involving strongly increasing coefficients. 
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Neglecting all w^s in the expansion of Wk($) = $ + W2 i T + w s^- • • • , we obtain for the 
expansions of the auxiliary functions f 1,2,3 as defined in Eq. fl28|) (recall that bi = \vi\V%&)' 



1=1 



(2i)\ 1 Zl ^(2i-l)\ 



;=i i=o ^ 

where i?2i denotes the Bernoulli numbers, and we define l/(— 1)! = 0. The r» are defined 
in Appendix |E| and are related to the group theoretical factors ^^ =1 _1 (i /2 ) 1 that occur in 
the expansions given above. Whereas the expansion of fa vanishes in the present approxi- 
mation, the expansion of its derivatives, as they occur in the last line of Eq. (pi]), must be 
retained: 

(d t - 4$<% + d)f 3 [ sw k , 



^•5^2^j_=5^^ a "" ,+1 ' (C - 2 ' 



tUj— >0 j = i 

D Cutoff functions 

In Eq. (|i~2f) , we introduce the cutoff function Rk(x) = xri-^), where r(y) is a dimen- 
sionless function of a dimensionless argument. For actual computations, we need the 
combinations h(y) and g(y) as well as their Laplace transforms h(s) and g(s) as defined in 
Eqs. (|nj) and (0). 

Instead of choosing a certain cutoff function by specifying r(y), we can specify a function 
h(y), or alternatively which fixes the remaining functions by virtue of Eq. (|TED; the 

direct connection between h(y) and g(y) can be formulated as 

vfg{v) = (g(v)-i)Ky). (D.i) 

A similar reasoning holds for a definition of the cutoff in Laplace space by specifying one 
of the functions h(s) or g(s), for which Eq. ( p.l[ ) translates into 



9(s) + s^-~g(s) = ~h(s) - [ dth(t)g(s-t). (D.2) 
as J 

These identities can be used to define a desired cutoff in its simplest representation without 
the need to specify the corresponding function r(y) explicitly; the latter might look very 
complicated. Of course, one has to take care of all the necessary conditions that a cutoff 



has to satisfy as listed in Eqs. (fL3|) and fll4|) . 



During the expansion of the propertime integrand in Sect. 0], we encounter the moments 



of h(s) and g(s) as defined in Eq. (Elf ). These moments can also be translated into a 
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momentum space calculation ( u y space"): 



oo oo 

f ds I /* 

h-j ■= I -^Hs) = — — J dyy^hiy), j > 0, 



o 

oo 



hj := [ ds s J h(s) = \im( -IV ( Hy), j>0 (D.3) 
J y-o \dyj 

o 

and equivalently for the g^s. 

In this work, the exponential cutoff is technically advantageous; all functions involved 
have a simple representation: 

1 y 



oo 



^0) = -E^ s - m )^' g(s) = S(s-1), (D.4) 

m=l 

where the s derivative acts on the remaining propertime integrand. For the moments 
required in d = 4, we find 

9j = !> 

h_ 2 = 2C(3) ~ 2.404..., (D.5) 
h j = B J> i = 1 > 2 , 
where Bj symbolizes the Bernoulli numbers. 



E SU(2) versus SU(3) 

Gauge group information enters the flow equation via the color traces. In Appendix |B|, 
we evaluated these traces formally by introducing the eigenvalue of (n c T c ) ab — > v\ y I = 
1, . . . , N 2 — 1. During the expansion of the right-hand side of the flow equation in Sect. f|, 
we encounter the following factors: 

JV 2 -1 

vf = n ai n a2 . . . n a ^ tr c [T {ai T a2 . . . T a2 % (E.l) 

i=i 

where the parentheses at the color indices denote symmetrization. For general gauge 
groups, these factors are not independent of the direction of n a . Contrary to this, the 
left-hand side of the flow equation is a function of \F^ V F^ U — > |_B 2 , which is independent 
of n a . Therefore, we do not need the complete factor of Eq. ( |F.1| ), but only that part of 
the symmetric invariant tensor trcfT 1 - 01 . . . T a2i ^] which is proportional to the trivial one: 

ti c \T^T a2 . . . T a2 ^} = n 5 {aia2 . . . 6 a2i _ ia2i) + . . . , (E.2) 
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where we omitted further nontrivial symmetric invariant tensors. These omitted terms do 
not contribute to the flow of Wk(9), but to the flow of other operators which do not belong 
to our truncation, e.g., operators involving contractions of the field strength tensor with 
the d a bc symbols. 

For SU(iV) gauge groups, we trivially deduce that 

r = N 2 - 1, n = N. (E.3) 

For the gauge group SU(2), all complications are absent, since there are no further sym- 
metric invariant tensors in Eq. (|E.2|) , implying 

7f= 2 = 2, z = l,2,... . (E.4) 

For the gauge group SU(3), we do not evaluate the Tj's from Eq. ( |E.2| ) directly; instead, we 
exploit the fact that the color unit vector can always be rotated into the Cartan subalgebra. 
For SU(3), we choose a color vector n a pointing into the 3 or 8 direction in color space, 
representing the two possible extremal cases: 

7V=3 _ 2 , 1 N=3 _ ^ (E z\ 

Note that their limiting behavior is rather different: for i — > oo, we find t^ =3 — > 2, but 

T i,S ^ U - 

The uncertainty introduced by the artificial n a dependence of the color traces is finally 
responsible for the uncertainty of our result for the SU(3) infrared fixed point. 



F Numerical computations 



Since the numerical evaluation of the anomalous dimension rj depending on the coupling 
G = g 2 /[2(Air) 2 ] as represented in Eqs. fl54|) and fl56|) is not straightforward, we mention here 
some details about the multidimensional integration and summation. We begin with the 
part rja. in Eq. (|54]): substituting S1/.S2 — * si, the $2 integral can be performed, resulting 
in the modified Bessel function K^ N 2_ 4 {2y/sl). Substituting furthermore t — > t/l, and 
defining the expressions 

: = £ ^TT^hl V K ^ = s '" 2/2_1 (F.i) 

1=1 



we obtain the representation 

oo 

64NG 



V. 



N=2 



T(W 2 -2)tt 2 



dtL(t) / ds 1 K(s 1 ) 



' f cGs.t 2 



1 (cGs x t 2 



(F.2) 
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Apart from an easily integrable singularity induced by L(t), the integrals are smooth 
and drop off exponentially for large t and si in the required G range. The sum S(q) defined 
in Eq. fl5"5|) converges quickly and an accuracy with error < 1% requires only (9(100) terms 
or less. The sum L(t) is rather slowly converging for small t, but the same accuracy can be 
obtained by including O(10 5 — 10 6 ) terms. Depending on the actual value of the arguments 
t and q, we adjust the included number of terms dynamically. 

For the part i] h , different complications occur. Beginning with Eq. (|56|) , we substitute 
s — ► styj cG / (2n 2 ) (and similarly for s') and find 



32NG 



T(3N 2 -2) 



Re 



Hit 



cG 
2^ 



;f.3) 



where we defined 



<*> = i/-(s) 

o 



1SS 



;f.4) 



The problem here is that the derivatives cannot be carried out numerically with a sufficient 
accuracy, but have to be computed analytically within the sum representation for S(-iss'). 
This implies that each term in the sum then consists of ~ 20 terms for SU(2) and ~ 50 for 
SU(3). This limits the generalization of the calculation to higher gauge groups for technical 
reasons. The remaining s and t integrations can easily be performed to a high accuracy. 
We estimate the total error of the numerical computation to be within a few percent. 
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